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1. Introduction 

I In the past ten years, optimization methods for image processing task have made a lot of 

I progress, thanks to the development of combinatorial methods (maximal flow/minimal cut, and 

other graph-based optimization methods — see for instance [10] , and [U [4] and the references 
therein). These methods are not new, the idea of representing Ising energies (i.e., discrete approx- 
imations of perimeters) on graphs and computing minimum points using maximal flows algorithms 
dates back at least to the 70s [13]. However, the evolution of computers and development of new 
algorithms [1], oriented towards specific applications, have contributed a lot towards the recent 
increase in activity in this field. In image processing, the idea is to regularize ill-posed inverse 
fyr") ' problems for finding sets (shapes) or partitions into labels of an image, by penalizing a discrete 

variant of their perimeter. We try to consider, in this paper, the most general energies which can 
be tackled by these methods, and even a little bit more: we consider discrete submodular energies, 
defined on discrete subsets of a finite lattice V C hZ^ , h > 0, for which it is known that polyno- 
mial algorithms do exist (see for instance [7] [TTl [T3] ) . We will show that, appropriately extended 
I into functions of general vectors in by means of the generalized coarea formula, these energies 

■ are, in fact, convex. This is already known (although our setting is a bit different, as well as our 

proofs which apply to other situations, including functionals defined in the continuous setting) in 
^> I discrete optimization, under the notion of Lovasz' extension [12j . 

^ • We will then study the continuous limit of our energies, as the discretization step h goes to (and 

the number of pixels/voxels in V to infinity), providing a very simple representation formula for 
the limit. In particular, it will be obvious from this formula that simple approximation procedures 
only provide "crystalline" energies, as already observed for instance in [3]. 

To be more specific, we consider in this paper an "interaction potential" F : {0, 1}^ [0, +oo), 
which is a nonnegative function of binary vectors of {0, 1}^, where S C is a finite (small) set 
of "neighbors" . We assume, in addition, that F satisfies the submodularity condition 

F{u Av) + F{u V v) < F{u) + F{v) 

for any u,v £ {0, 1}^, where {u A v)i := mm{ui,Vi} and {u V v)i := ma.x{ui,Vi}. Defining, for 
X G and u a real-valued function, the vector u[x + hY,] — {u{x + hi)i^s) G M^, we will study 
the asymptotic behavior as ft. ^ of functionals of the type 

(1.1) ME,n):^h''-' F{xE[x + hE]), 

xe/''(n)n/iZ" 

where here, £7 is a bounded open subset of and is a discrete subset of the discrete lattice 
V = hZ^ nil {E is also identified to the union of the cubes = a; -I- [0, , x £ E, and xe is 
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its characteristic function). The notation stands for the points x such that x + hH C il, so 

that the sum in p.ip involves only the nodes x G /iZ^ such that x + hi] C fl. 

The functional (|l.ip is a sort of nonlocal anisotropic discrete perimeter of E. In fact it penalizes 
the boundary of E in a nonlocal way, since an interface at the boundary with a vertex x interacts 
with the behavior of E on the cubes with vertices x + hT,. The nonlocality vanishes as ft, ^ since 
its radius of action is given by h diam{T,). The anisotropy is introduced by the function F, which 
can weight interfaces with various orientations in different ways. 

The main result of the paper concerns the asymptotic behavior of the discrete perimeters (jl.ip 
as /i in the variational sense of F-convergence (see Section [5] for the definition) with respect to 
the L^-topology on the family of discrete sets (that is convergence of characteristic functions). 
Under mild assumptions on F and fl, we prove that (see Theorem 14. 2p the discrete perimeters 
F-converge to the continuous anisotropic perimeter which for a sufficiently regular set E (a set 
with Lipschitz boundary for instance) is given by 

(1.2) J{E,n)^[ F{uE-T.)dA, 

JdE 

where ve is the inner normal at the boundary and {ve ■ S) = {ve ■ y)y£^- More precisely the 
class on which the F-limit is defined is given by the family of sets with finite perimeter in fi. As 
a consequence, for a general set E, the boundary involved in the functional (|1.2p is the reduced 
boundary d*E, the inner normal ve is intended in a measure theoretical sense (see Section[2|), and 
the area measure dA has to be replaced by the {N — l)-dimensional Hausdorff measure TC^"^. 
The function F appearing in (|1.2p is the extension to of the submodular function F by means 
of the formula 

/+00 
FiX{u>s})ds, 
-oo 

where {u > s} := {a; G E : u{x) > s}. Formula (jl.3p is a coarea formula for the function F since it 
relates the value F{u) to the behavior of F on the "boundary" of the level sets {u > s} (compare 
with equation (|2.3p which gives the classical coarea formula for functions of bounded variation). 

In view of the result on dicrete perimeters, we obtain a F-convcrgence result for the functionals 
J/i(-,r2) extended to the class of piecewise constant functions u relative to the grid hZ^ . More 
precisely we consider u of the form 

u = ^ a^XQ^, G R 

and 

(1.4) Jh{u,n)^h^-^ F{u[x + hl]). 

a;G/''(0)nhZ" 

As the functional (jl.ip could be thought as a discrete perimeter, the functional (|1.4p could be 
considered as a sort of discrete total variation of the function u. Clearly it inherits the nonlocal 
and anisotropic features of the discrete perimeter. We show that the F-limit in the L^-topology is 
given by the anisotropic total variation 

(1.5) j(^u,n)= I^f(^^^-I^ d\Dul 

where u belongs to the space BV{fl) of functions of bounded variation (see Section[2]). Du/\Du\ G 
gAf-i c;igjioting the Radon-Nikodym derivative of Du with respect to its total variation \Du\. 
This F-convergence result is a simple consequence of the result on sets and of the fact that the 
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functionals satisfy the generalized coarea formula 



/>+oo 



(1.6) 



Jh{u,Q)= Jh{x{u>s},^)ds, 



so that the behavior of Jh on piecewise constant functions is completely determined by the discrete 
perimeters for sets. We infer the result from general properties of functionals on L^{n) that 
satisfy a coarea formula like (|1.6p . which we study in Section [3l This class of functionals. denoted 
by GC(r2), was investigated by Visintin [151 116| in connection with phase transition problems. 
As a consequence of our F-convergence result, the discrete total variations p.4|l can be used 
to approximate Total Variation Minimization procedures in image denoising involving (|1.5p (see 
Corollary [43]). 

The paper is organized as follows. Section [2] contains the notation employed in the paper, 
and some basic facts concerning sets with finite perimeters, functions of bounded variation and 
F-convergence. In Section[3]we consider the class GC{^1) of functionals on L^{n) which satisfy the 
generalized coarea formula: In particular we prove that GC(f2) is closed under F-convergence, and 
that the limit can be recovered by the behavior on characteristic functions of Borel sets. Section 
E] contains the main F-convergence result formulated for the discrete total variations p.4p . We 
exploit the reduction to the class of discrete sets, and Subsections 14.11 and 14.21 contain the proof 
of the two inequalities characterizing F-convergence for study the discrete perimeters. 



Let A be an open subset of M^. We will say that A has a continuous boundary if OA can be 
covered by finitely many balls B such that, in each ball, B O A is the subgraph of a continuous 
function (after an appropriate change of coordinates). If these functions are Lipschitz continuous, 
we say that A has a Lipschitz boundary. 

For any p Cz [1, +oo[ we will denote by LP{A) the usual space of all p-summable functions on 
A, and by L^{A) the space of measurable functions on A which are essentially bounded. Given 
u, V e L^{A), we set 

(2.1) M A V := min{u, w} and w V w := maxjii, w}. 

In the following, we recall some basic facts concerning function of bounded variation and sets 
with finite perimeter which we need in the following sections, together with some basic definitions 
and results concerning F-convergence. 

Functions of bounded variation and sets with finite perimeter. For an exhaustive treat- 
ment of the subject, we refer the reader to 

We say that u has bounded variation in A and we write u g BV{A) if u G L^{A) and 



|Dtt|(>l) is referred to as the total variation of u. 

li E C A is a Borel set, we say that E has finite perimeter in A if xe £ BV{A), and we set 



2. Notation and preliminaries 



(2.2) 




Per{E,A) ■.= \Dxe\{A). 



Per{E, A) is called the perimeter of E in A. It turns out that 



Per{E,A) =n^-\d*EnA), and Dxe 



veH 



\_d*E. 
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where d* E denotes the reduced boundary of E, which, up to a Ti.^ ^-neghgible set, coincides with 
the (larger) set of points x such that there exists a unit vector with 

E — X 

— > {y e : ?/ • i^Eix) > 0} as £» ^ in LI^{R^). 

The unitary vector ve{x) is usuaUy referred to as the interior normal to E at x. Ti^~^ denotes 
the {N ~ l)-dimensional Hausdorff measure, which is a generahzation to arbitrary sets of the usual 
[N — l)-area measure. The points of d* E are also called regular points of dE. 
If 7i G BV{A), the following coarea formula holds: 

/+00 ;'+oo 
Per{{x&A:uix)>s},A)ds= / \Dx{u>s}\{A) ds. 
-oo <y —oo 

Finally we recall the following compactness result (which is a variant of Rellich's theorem). If A is 
bounded and with Lipschitz boundary, and (w„)„gN is a sequence in BV{A) such that + 
|I?u„|(A) is bounded, then there exist a subsequence (u„j.)fegN and a function u e BV{A) such 
that 

w„fc ^ u in L^{A) 

and 

\Du\{A) < liminf |i:>w„J(A). 

k — 'OO 

F-convergence. Let us recall the definition and some basic properties of De Giorgi's T -convergence 
in metric spaces. We refer the reader to [HIIS] for an exhaustive treatment of this subject. Let (X, d) 
be a metric space. We say that a sequence F„ : X [— oo, +oo] F-convergcs to F : X ^ co, +oo\ 
(as n — > oo) if for all u £ X we have 

(i) {T-lim inf inequality) for every sequence (w„)„gN converging to u in X, 

liminf F„(m„) > F{u); 

n — >OD 

(ii) {T-lim sup inequality) there exists a sequence (w„)„gN converging to u in X, such that 

limsupF„(M„) < F{u). 

n — 'OO 

The function F is called the F-limit of (i^„)„gN (with respect to d). Given a family (i^/i)h.>o of 
functionals on X , we say that F^ F-converges to as /i ^ if for every sequence hn ^ we have 
that Fh„ F-converges to F as n ^ oo. 

F-convergence is a convergence of variational type as explained in the following Proposition. 

Proposition 2.1. Assume that the family {Fh)h>o T-converges to F and that there exists a 
compact set K C X such that for all h > 

inf Fhiu) = inf Fhiu). 

Then F admits a minimum on X , mix s- minx F as h ^ 0, and any limit point of any 
sequence {uh)h>a such that 



\im(Fhiuh)- inf Fh{u)) =0 



is a minimizer of F. 
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3. Generalized coarea formula 
In the following, let il C M.^ be an open and bounded set. 

Definition 3.1. Let J: L^{n) [0, +oo] be a proper functional. We say that J satisfies the 
generalized coarea formula if for every u € L^(f2) 

/+00 
JiX{u>s})ds, 
-oo 

with the convention J{u) = +oo if the map s t-^ Jixiuys}) *s not measurable. We denote by 
GC{fl) the class of functionals satisfying (|3.ip . 

The class GC{rt) has been introduced by Visintin [15] and investigated, in the discrete case, 
by ChamboUc and Darbon [5] . In a slightly different setting, the formula (|3.ip is a variant of the 
extension introduced by Lovasz in [12| and well-known in combinatorial and linear optimization. 

An example of functional satisfying p.ip is given by the total variation (|2.2p in view of the 
coarea formula (|2.3p . Other examples are treated in |16j : 

J(m)=/ \u{x)-u{y)\\x-y\-^^+'^Uxdy, Vr G (0, 1) 
JnxQ 

and 

J{u) = / (ess sup u — ess inf dxdh, Vr G (0, 1). 

Jnxm+ Bhix)nn Bh{x)nn 

The next Proposition contains some elementary consequences of formula (|3.ip . 

Proposition 3.2. Let J G GC{fl). Then for every u G L^{Vl) the following facts hold: 

(i) J{Xu) = XJ{u) for every X > 0; 

(ii) J{u + c) = J{u) for every c G K; 

(iii) J(c) = for every c G M. 

Moreover if J is convex, for every u, v G L^{n) we have 

(iv) J{u Av) + J{u y v) < J{u) + J{v). 

Proof. Let A > 0, M G L'^{n) and c G K. Then 

/ + 00 /•+00 
J{X{u>j.))d'S ^ ^ J J{x{u>t})dt ^ XJ{u) 

and 

/+00 ;'+oo 
J{X{u>s^c})ds^ j J{X{u>t})dt = J{u) 
-oo J ~oc 

so that (i) and (ii) follow. 

Let us prove (iii). In view of (ii), it suffices to show that J(0) = 0. Suppose by contradiction 
that J(0) > 0. Then for every u G i°°(rj) we have 

/-(-oo r+oc 
J{X{u>s})ds> / J(0)ds = +oo. 

-oo CSS sup u 

n 

By the generalized coarea formula, we deduce that J{u) = +oo for every u G L^(il). But this is 
against the fact that J is proper, so that point (iii) is proved. 
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Let US show (iv). Since we have 



J{X{u>s} A X{v>s}) + J{X{u>s} V X{v>s}) = / •^(X{X{u>.}AX{„>.}+X{„>.}VX{„>.)>t}) 







+00 



J(X{xi }AX{„>,}+X{„>s}VX{„>s}>t}'"^''' 



by the generahzed coarea formula (|3.ip we get 

'^(X{tt>s} A X{v>s}) + J{X{u>s} V X{v>s}) = J{X{u>s} A X{v>s} + X{u>s} V X{v>s}) 

= J{X{u.>s} + X{i»s})- 

Notice that if J is convex, by point (i) we deduce that J is subadditive. Then we obtain 

(3.3) J{X{u>s} A X{v>s}) + J{X{u>s} V X{v>s}) < J{X{u>s}) + J{X{v>s})- 

Observe that for any s S R we have {u A v > s} ^ {u > s} O {v > s} and {uV v > s} ^ {u > 
s}U {v > s} so that 

X{uAv>s} = X{u>s} A X{i>>s} and X{uwv>s} = X{u>s} V X{v>s}- 
We conclude by (|3.3p 



f* + QO 

J{uAv) + J{u\/v)= I [J{X{uAv>s}) + J{X{uwv>s})]ds 

-00 
+00 

[J{X{u>s} A X{d>s}) + J{X{u>s} V X{f>s})] rfs 

— 00 

-\-oo />-\-oo 

JiX{u>s})ds + / JiX{v>s})ds = J{u) + J{v) 



< 



so that (iv) follows and the proof is complete. □ 

We will need the following Lemma concerning the approximation of Lebesgue integral by means 
of Riemann sums. 

Lemma 3.3. Let f e i^(M), t e]0, 1[ and let us set 

T7 ^ ^ 



n 



Then up to a subsequence we have 

/+00 
f{T)dT for a.e. t e]0, 1[. 

Proof. For any 77, G N we easily get 
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Then for t elO, 1[ we have 



+00 



/(T)dr-s„(t) 



dt 



+ 00 



n ^-^ \ n I 

kez ^ ' 



But 



1 /.I 



JO 



fc + r 



71 ^ 

feGZ 



n 

1 (-1 r 



Jo 



/ 



-/ 
k + r 



k + t 



f 



k + r 
n 

< ■ 

dr dt 

k + r + t 
n 

t 
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E 



1 rl 



k + t 



dt 



dt 



JO 



n 

k + r 
11 



k + t 



dr dt. 



f 



k + r + t 



f 



k + t 



dr dt 



dt 



LH 



/(• 



/ • + - 

n 



dr. 



/(•)-/ • + - 
\ n 

The last terms tend to zero by continuity of the translation operator in L^(M). We conclude that 

f +00 

f{T)dT inii(0,l) 

so that, up to a subsequence, pointwisc almost everywhere convergence follows. □ 

In view of (|3.ip . functionals in the class GC{U,) are completely determined by their behavior on 
characteristic functions of Borel sets contained in fi. The next result gives a sufficient condition 
for the convexity of lower semicontinuous functionals in GC(f2) in terms of the submodularity 
property (iv) of the previous Proposition only on characteristic functions. 

Proposition 3.4. Let J G GC{^) he a lower semicontinuous functional such that 

(3.4) JixEnE' ) + JiXEuE' ) < J{xe) + J{XE' ) 
for every pair of Borel sets E, E' in il. Then J is convex. 

Proof. Since by Proposition 13. 21 J is positively one- homogeneous , it is sufficient to show that 

(3.5) J{u + v) < J{u) + Jiv) 
for any m, u G L^{fl). 

We claim that the following representation formula holds for every function u which is positive, 
bounded and with integer values: 



(3.6) 



J (u) = min I ^ J (xEi) ■■m>l,u^^XEi 



In order to prove p.Sp . we can clearly assume that J{u) < +00 and J{v) < +00. Hence by p.ip 
the maps s 1— *■ J(x{u>s}) ^^'^ J{X{v>s}) belong to L-'^(R). 

Firstly let us assume < u < 1 and < v < 1. For every n G N, n > 0, let us set for t g]0, 1[ 



-In := - E 



feeN 



and 



fcGN 
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By applying Lemma 13.31 wc can choose t g]0, 1[ in such a way that 

1 



fceN 

and 



1 /"^ 



fcef 

By construction u„ u and ^ in L^(fi). Then by positive homogeneity, and assuming the 
representation formula (|3.6[) holds, we get 

\ 1 NT- / \ / 

} 



J{Un+Vn) = J I^^X{u>!s±±}+X{y>h±lj J < ^^Yj'^{^u>^}) +j{x{y>h±l 

\ fceN " " / fceN 



The right-hand side converges by construction to J{u) + J{v), and thus, by lower semicontinuity 
of J, we have that (|3.5p follows. 

In the case m < u < M and m < v < M, one can easily show again that (|3.5|1 holds by 
considering the functions {u — m)/{M — m) and {v — m)/{M — m), and taking into account the 
general properties of J. 

Finally, for m, t> € L^{^) and for T > 0, let us consider ut := -T\J u A T and wt := -TV v A T. 
Since mt — > u and iij^ — > w in L^(f2), by the lower semicontinuity of J we obtain 

J{u + v) < liminf J(ut + wt) < liminf( J(ut) + J(ut)) < lim sup J(ut) + lim sup J(wt) 

T-»+oo T-»+oo T^+cx) T^+oo 



= ^1™ / "'^(X{u>s})c^s + ^lim / •/(X{i.>s})ds = -/(w) + J(w), 

1 > -|- OO J rp 1 > -\- OO J rp 

SO that p.Sp follows. 

In order to conclude the proof, we have to check claim (|3.6p . Let 7\f := maxu. Since u is 
positive and integer valued, we can write u = X{u>i}- For any i £ {1, . . . , M} we have 

J{X{u>s})ds = J(x{«>i}) 

i-i 
so that 

J{u)= i J{X{u>s})ds = ^J{x{u>^})■ 
''^ i=l 

Hence 

{m m ^ 

In order to prove the opposite inequality let u = X-E; for some Borel set Ei C_ fl and m > 1. 

Observe that for any s G {1, . . . , m} with r ^ s wc also have 



From (|3.4p we get 



CONTINUOUS LIMITS OF DISCRETE PERIMETERS 



9 



Then by induction it is easy to see that 

f m m ^ 

inf < ^J{XE,) ■■m>l,u = ^XEi\ 



M 



> inf <^ ^ J{XE.) ■■m>l,u^Y.^E,^EiDE2D---DErn>^J2 JiX{u>^}) = J{u). 

I 1=1 1=1 J i=l 

Hence claim p. 61) holds true, so that the proof is concluded. □ 

The following Proposition deals with the stability of the class GC{Q) with respect to the F- 
convergence. 

Proposition 3.5. Let (Jn)neN be a sequence of convex Junctionals in GC{^) such that there 
exists a functional J defined on characteristic functions of Borel sets which satisfies the following 
conditions: 

(a) for every Borel set E Q and for every sequence of Borel sets (£'„)„gN contained in Q 
such that xe„ Xe in L^{p) we have 

J{xe) < liminf J„(x£;„); 

n — 'OO 

(b) for every Borel set E C Q there exists a sequence of Borel sets (£'„)„gN contained in Q 
with xe„ —* Xe in L^iVl) and such that 

limsup J„(x_E„) < J{xe)- 

n — >oo 

Then setting 

(•-f oo 



J(u) := / J{X{u>s})ds, 

J —oo 

we have J € GC{Q) and the sequence {Jn)n&i T -converges to J in the -topology. 

Conversely let {Jn)n£n be a sequence of functionals in GC(fl) which T-converges to a proper 
functional J: L^(fl) [0, +oo]. Then J G GC(f2) and its restriction J to the family of charac- 
teristic functions of Borel subsets of satisfies conditions (a) and (b). 

Remark 3.6. Notice that it follows that for convex functionals in GG{fl), the F-convergence 
is equivalent to the F-convergence on the corresponding (submodular) set functions, that is, the 
restriction to characteristic functions of the original functionals. However, the last statement in 
Proposition 13 . 51 is also true without assuming any convexity of the functions J„. Notice that there 
exist functionals in GG{Q) which are lower semicontinuous but not convex, so that convexity 
cannot be gained by relaxation. (It suffices to consider functionals of the form (|4.4p with fl and 
h chosen in such a way that the summation involves only one square, and the function F is not 
submodular on binary vectors.) 



Proof of Provosition \3.5l Notice that J is, by construction, lower semicontinuous on characteristic 
functions, so that the map s t—f J{X{u>s}) is measurable for every u G L^{fl). Hence the definition 
of J is well posed. 

In order to prove the F-convergence result, we need to check F-lim inf and F-lim sup inequalities 
(see Section [2]). Let us start with the F-lim inf inequality. Let u„ ^ u in i^(fi). Up to a 
subsequence, we can assume that X{u„>s} ~* X{u>s} in L^{^) for a.e. s G K. By Fatou's Lemma, 
the generalized coarea formula (j3.ip and assumption (a) we get 

/+00 I'+OO 
liminf J„(x{u„>s}) ds > JiX{u>s})ds = J{u) 

-oo J-oo 

so that the F-lim inf inequality follows. 
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Let US come to the F-lim sup inequality. Wc can clearly assume that the map s i— > J{x{u>s}) 
belongs to L^(R). Notice that the subspace given by (finite) linear combinations of characteristic 
functions is dense with respect to the energy J. In fact, if < w < 1, by Lemma 13.31 we can choose 
t e]0,l[ such that 

keti 

and (since J(0) = J(l) = 0) 

J^^I^lom^'^ V^b>^}) ^ %{n>.})o's=/ J{X{u>s})ds. 

The case m < u < M follows considering the function {u — m)/{M — m). Finally, for u G L^{n), 
let us set ut —T \/ u A T ^ u for every T > 0. Since ut — > u in L^(Q) and J{ut) —> </(«) 
as T ^ +00, the density in energy follows by a diagonal argument. By general results of F- 
convergence, it suffices to prove the F-lim sup inequality for u equal to a linear combination of 
characteristic functions. Since J is invariant under translation, it is not restrictive to assume that 
^ ~ SfeLi o^kXEk with flfc > for every k = 1, . . . , m, and Em Q E„i-i ^ ■ ■ ■ Q Ei. In this way 
we have 

By condition (b), we can find Borel sets EJ^ such that 
as n — > oo, and 

limsup J„(x£;") < J(xEk)- 



n — >-oo 



Setting Un X^^Li^fcX-E^? have Un u in L^{Q). Since Jn is convex and positively one- 
homogeneous, and hence subadditive, we deduce 

m m 

limsup J„(m„) < ^Ofclimsup J„(x_E^) < ^a/cJ(XBj = J{u), 

k—1 k=l 

SO that the F-lim sup inequality is proved. 

Finally, the fact that J S GC(ri) follows since J and J coincide on characteristic functions. 
The proof of the first part of the Proposition is thus complete. 

Let us come to the second part. Clearly J satisfies condition (a). In order to prove condition 
(b), let i? be a Borel subset of ft, and let u„ S L^{n) be such that w„ xe in L^{^1) and 
limsup„^3o Jn(un) < J(xe)- Since for any 5 € (0, 1) 

^1-5 

Jn{Un)> J JniX{u„>s})ds, 

there exists s„ G ((5, 1 — S) such that 

(3-7) JniX{u„>s„}) < 1 ^2(5 

Let us set :— X{u„>s„}- We have clearly that Xe'' ^ Xe in L^{^) and by p.7p we deduce 

limsup J„(x£;i) < jlimsup J„('u„) < ^ ^r- JiXE)- 

Let US choose now d = 1/m. There exists Um such that for every n > rim we have 

llXiji/- -XeWl^ < l/m 
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and ^ 

MXpU^) < Z TTT-JiX 



E 



t+1 



Moreover we may assume that '\ oo. If we set '■= E}J^^ for < n < rim+i, we have that 
(i?„)„gN is the recovering sequence for which the r-Um sup inequahty holds. 

Finally, the fact that J G GC(fi) follows now from the first part of the Proposition, and this 
concludes the proof. □ 

4. Discrete approximation of anisotropic total variation 

Let > 1 and S C be a finite set, and let F : {0, 1}^^ — > [0, +oo) be a nonnegative 
submodular function, i.e. F{u /\ v) + F{u V f ) < F{u) + F{v) for any u,v E {0,1}^, with 
F(0) = F{xs) ~ 0. We extend F to all vectors u E into a convex function by letting (see 
Proposition 13. 4p 

/ + 00 
Fix{u>s})ds 
'OO 

where {w > s} := {x e S : u{x) > s}. We let 

(4.2) := max|x| 

and we assume, in addition, the following coercivity assumption: 

(A) S contains and the canonical basis {ei)fLi of M^, and there exists c > such that for 
any u g M^, 

N 

F{u) > \u{e,) - u{0)\. 

i=l 

Notice that (|4.ip is a discrete version of the generalized coarea formula (|3.ip . 
Given /i > and a; G /iZ^ let us set 

(4.3) :=a; + /i[0,l[^. 
Let Vh denote the space of functions u : R such that 



u 

Notice that we have Vu C i;i„^(R^). 

Let C be an open and bounded set. We denote by Vh{^) the restriction to of the 
functions in Vh- Let moreover denote the set of x € such that x + ZiE C fJ. We consider 

the functional J/j(-,fi): i^(ri) [0, +oo[ defined as 

f/i^-i F{u[x + h^) ifueVhin) 

(4.4) J/i(ii, fJ) := < 2;G7''(n)n/iZ" 

[+00 if u e Li(f7) \ n(17), 

where for any x G /''(fi), 'u[a; + /iS] is the vector {u{x + hy)y^s) of R^. 

The aim of this section is to study the asymptotic behavior of the functionals J/i(-,i7) as the 
size mesh h vanishes: it is expected that they approximate some anisotropic total variation. The 
following Proposition shows that the functionals Jh{-,^) satisfy the generalized coarea formula 

m- 

Proposition 4.1. The Junctional J h{- is convex and belongs to GC{n). Moreover, there exist 
C2 > Ci > such that for any open sets A, B with A CC B CC fJ, and for any u £ Vh{^), we 
have, if h is small enough, 

(4.5) Ci\Du\{A) < Jh{u,B) < C2\Du\{n). 
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Proof. From (|4.ip . wc get that also Jh satisfies p.ip . The submodularity of F yields (|3.4p . hence 
J/j is convex. 

To show the estimate (|4.5p . it is enough to assume that u £ V/j(ri) is a characteristic function 
(the general case then follows from the coarea formula). In this case, the left hand side inequality 
follows from assumption (A), while the other follows from the fact if that F{u[x + /iS]) > for 
some X G I^{B) n ft.Z^, then u takes different values on the set x + /lE so that its variation on 
B{x,p-s:h) (where ps is given by (|4.2p ') is at least h^~^. □ 

For every i/ E we set 

F(i..E) :=F((t..y),es) 
The main result of the paper is the following. 

Theorem 4.2. Let n C be a bounded, open and Lipschitz domain, and let Jh J/i(-,f7) be 
defined as in (|4.4[) for h > 0. Then the family {Jh)h>o T-converges in the -topology as h —^ 
to the functional J : L^(Jl) [0, +oo] given by 

Du 




, , . , , ,1] d\Du\ ifue BV(n) 
(4.6) J{u,n)={Jn \\Du\ '11^ ^ ' 



ifue L\n) \BV{n), 

where for u G BV{^) the function Du/\Du\ stands for the Radon-Nikodym derivative of Du with 
respect to its total variation \Du\. 

Since J(-,ri) satisfies the generalized coarea formula (see Proposition 13. 5p . we can also write 
for u g BV{n) 

/ + 00 
Pers,F({M > s},Q.)ds, 

where for any finite-perimeter set E 'm ^ 

(4.8) Per^,FiE,n) = [ F{ve ■ T.) dU^'^ . 
In particular 

J{XEM^Per^AE,n) 

for any finite-perimeter set E in 57. 

The following Corollary is a consequence of the F-convergencc result of the previous Theorem. 

Corollary 4.3. Let Vl C be a bounded open set with Lipschitz boundary, g € L°°{^), and let 
Uh be the solution of 

(4.9) min Jh{u, 17) + ||u - .g|||2(m. 

Then Uh converges in L^(f2) for h —> to the minimizer u G BV{fl) of 

(4.10) min J(u,r!) + ||u-g|li2(a), 
where J is the V -limit of the family { Jh)h>o given by (j4.6p . 

Proof. Without loss of generality we can suppose Uh S Vh{fl) for any h > 0. Moreover, as a 
consequence of the submodularity property (iv) in ProDOsition l3.2[ we have that the functional Jh 
decreases by truncation. This entails that ||uh||oo < llslioo for every h > 0. 

Taking into account (|4.5p we get that the total variation of Uh is uniformly bounded. By 
compactness in BV, we deduce that there exists u e BV{il) and a sequence h^ ^ such that 
Uhk —> u in L^(n). The convergence is indeed in every L^ for every 1 < p < +oo since (u/i)h>o is 
bounded in L°°{n). 
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The fact that the hmit u is a minimizcr of (|4.10p is a consequence of Proposition 12.11 Since 
this minimizer is, in fact, unique, we conclude that the entire family (it/i)/i>o converges to u as 
h^O. □ 

Remark 4.4. Notice that equality (|4.7p implies that the F-limit J{u, J7) satisfies (|4.6p for every 
u € BV{il). In fact a direct computation shows that (|4.6p holds for simple functions. The 
extension to the whole BV{fl) follows by a density argument. Let u E BV{U), and let (m„)„(=n 
be a sequence of simple functions converging in L^(ri) to u and such that \Dun\{il,) \Du\{n) as 
n — > oo. From Rcshctnyak continuity theorem (see [2] Thm 2.39]) we deduce that 

i^4/{w:i ■ "i^-i = /„ <S ■ "i^-i- 

Since iZJunKil) jZJuKri), by coarca formula in BV{il) wc get, up to a subsequence, Per{{un > 
s}, H.) Per{{u > s}, il) for a.e. s G R. Using again Reshetnyak continuity theorem (applied to 
the measures v dH^^^\—d*{un > s}), we deduce that 

lim Pers f({u« > s}, O) = lim / Fii^ ■ S) dH^-^ 

' Ja*{«„>s}nn 

= / F{u- E) dn^-^ = Perj:,F{{u > s}, n). 

Jd*{u>s}nn 

By the generalized coarea formula for J and the Dominated Convergence Theorem wc conclude 
that 

lim J(u„,r2) — J(m, f2) 

n — *oo 

so that in view of (|4.11[) . the representation (|4.6p is proved. 



The rest of the Section is devoted to the proof of Theorem 14.21 In view of Proposition 13.51 and 
of Remark 14.41 in order to study the F-limit of the family {Jh)h>o we can consider the restriction 
of Jfi to characteristic functions of sets, and show that it F-converges to the anisotropic perimeter 
given by 

By definition, we need to show that given any sequence hm i 0, we have for any Borel set 

• if x_B„ e 14^ (fi) converges to xe in L^{^), then 

(4.12) liminf A,„(XB„,^^) > J{XE,n); 

m — »oo 

• there exists a sequence (£'m)meN with XEm ^ ^h^i^) such that xe^ ~^ Xe in L^{^) and 

(4.13) limsupA„(XB„,l^) < J{XE,^). 

m — 'OO 

We prove inequalities (|4.12p and (|4.13p in subsections 14.11 and 14.21 respectively. We will use the 
following "continuous" variant of J^, defined on any function and not just on piecewise constant 
functions of the class Vh'- we let, for any u E L^{n), 

(4.14) Jfi{u,n) -.^ ^ [ F{u[x + h^)dx. 

Let Qi, is the open unit cube centered in with a face orthogonal to v, and 

(4.15) := {xeR^ ■.x-i^> 0}. 
We have the following result: 

Lemma 4.5. There holds 

(4.16) F{i.-E) = lim j;;(x/„,Q.) 
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Proof. One has 

Now, letting u{x) := v ■ x, we have, for any s G M and x £ , xi„ [x + /iS] = X{u>s} [x + + hJ^], 
so that we may write 



1/2 



1 



-1/2 \ h J Ih(Q^ 



(4.17) Jt{xi..Qu) = / J r L,_ ,-P^(X{u>.}[a^ + '^^ + /»S])c«a:J 

,1/2 



/ (1/ F(x{„>4[2/ + /iS])d2/) ds. 

J-l/2 V" / 



Now, as soon as \s\ < 1/2 — p^h, 



I F{x{u>s}[y + hY.])dy = / F{x{u>8}[y + hY.]) dy , 

and it follows from (j4.17p and the co-area formula that 
(4.18) JfXxi^.Q.) - Jtiu^Qu) + eu 

where the error > is 
1 



h 



5-ps/l.<|s|< 



f / F{x{u>s}[y + hY])dy- / F{x{u>s}[y + hY])dy \ ds. 

5 \"'-f''(Q.)+^i' Ji'^iQu) / 



One easily checks that et < 2(maxx; F)p^h — > as /i ^ 0, and (|4.16p follows from (|4.18p and the 
observation that 

r,{u,Q.) ^ F(z..I])^^^^ ^ F(^.S) 

as h^O. □ 

4.1. r-liminf inequality. We prove here the inequality ()4.12p . First of all, we may assume 
(upon extracting a subsequence) that the liminf in (|4.12p is a limit, and, also, that it is finite 
(otherwise there is nothing to prove), so that in particular sup,„ J/j^ (x_b„ , i^) = C < +oo. For 
any A CC il, it follows from (j4.5p that for any m large enough, Ci|Z)x£;„ l(^) < C so that (by 
lower semicontinuity of the total variation) also E must have finite perimeter in fl. 
We consider the non-negative measures 

p„i = h'^^^ ^ F{xEAx + hra^)5x, 
xe/'"n(n)n/i,„z" 

such that JhmiXEmT^) — /^rn(^)- Since they arc uniformly bounded, we may also assume that 
there exists a measure p such that /i™ ^ /i as measures. We therefore have 

(4.19) p{n) < limini JhJxE,^,^) , 

m — >oc 

hence the thesis follows if we show that p > F{ve ■ J^)TC^~^\—d* E as measures. 

It is therefore enough to compute the Radon-Nikodym derivative of the measure p with respect 
to H^^^l—d^E, and to show it is above F{i'e ■ S). We know from [51 Thm. 5.52] that it is given, 
for Ti^'i-a.e. x e d*E, by 

+ rQ^) 

rZn^^^{x + rQ,)r]d*E) 
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for any v G with unit norm, where Qy is as before the open unit cube centered in with a 
face orthogonal to v. In particular, at a regular point .tq of d*E we can choose v = v^ix^) (the 
inner normal to at xo) and the limit becomes 



(4.20) 



i = hm 



r^O r 



JV-1 



Let us now show that I > F{v ■ E). We can assume ^ < +oo. Notice that since xq is regular, we 
also have 

/ \xiAy-xo)-XE{y)\dy ^ o. 



lim — TT 



r— »o r 



N 



For a.e. r > (small), we have 

^(xo+rQ^) = lim ^lm{xo+rQ^) , 

ni — »oo 

and 



WiAv - xo) - XE{y)\dy = lim 



xo+2rC, 



WiAy-xo) -XE„Ay)\dy- 



xo+2rQ^ 



Hence, using a diagonal argument, there exist sequences m„ and r„ such that hm^/vn — * 0, 

A*rn„(a;o + rnQ^) 



lim 

n — >QC 



and 



lim 



1 



\xiAy~XQ)-XErn„{y)\dy = o. 

For each n, we now make the change of variable y ^ xq + rnZ, and we define E'^ — {E„ 
xo)/rn C (fl — xo)/rn. It follows 

N-1 



(4.21) 

and 
(4.22) 



= lim 



zeQ„n[(/i„„/r„)Z"-2-o/r„] 



lim 



IXiA^) ~XE'^iz)\dz = 0. 



2Q„ 



We let h'j^ = hra„/rn (which goes to 0), and let 0„ = xo//im„ — [xo/hmr,] (the vector whose 
coordinates are each the fractional part of the corresponding coordinate of a;o//im„)- The limit 
in (|4.2ip becomes: 

E F{xE'^z+Kn). 

zeQ„n/i'„(z«-e„) 
Letting E[[ := S'^ + we clearly still have 



(4.23) 



lim h': 

n — >oo 



(4.24) 



lim 



\xi,Az) - XE''{z)\dz = 0. 



From now on, to simplify, we will denote E[[ and h'^ by, respectively. En and ft,„. We consider 
a basis {vi^ . . . ,vn) of such that each i/j is orthogonal to a face of Q,j, and with i^w = v. We 
choose 770 > (small). Writing x = sj^i + a;' = (s,^:') with x' • i^i = 0, we have (using Fubini's 
Theorem) 



(4.25) 



1 



1/2-r, 



Jl/2-rj-psh„ J [-1/2,1/2]" 



= PS 



Ix/. (2;) - XE„ ix)\ dx' ds 

1/2 n 

1/2-i-iQ-pT.h^ J [-1/2, 1/2]'' 



drj 



\xiAx) - XeSx)\ dx'ds 



which, thanks to ()4.24p . goes to 0. 
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Hence, up to a subsequence, we may assume that for a.c. rj G [0,770], 



1 



1/2-7, 



(4.26) — / / \xiAs^i+x')-XEAsi'i+x')\dx'ds ^ 

as n — > 00. The same holds if we replace in (j4.26p vi with any of the j/^, i ~ 2, . . . , — 1, or if we 
replace the interval [1/2 — ry — ps^-n, 1/2 — r/] of integration in the first integral with the interval 
[—1/2 + 77,-1/2 + 77 + pshn]- Let us therefore choose a 77 > such that all the above mentioned 
limits are 0, and let QJ^ = (1 — 27])Q^. 

We now extend periodically En n Q2 in the directions z^i, . . . ,J^Ar_i, into a set En C i^-^ + 
[—1/2 + 77, 1/2 — -qlv, by letting x £ En if and only if a; = a; + (1 — 2rj) X^ila^ ^i^i f^'' some point 
X e En C] Q]!, and integers ki S Z. It follows from (|4.23p that 

(4.27) i > liminf J,^Jxb ,Q2), 
where Jjl is as in f|4.14p . and (recalling (|4.5p ). 

(4.28) suplT^x^ \(Q2) < snp\DxEj{int{Q"l)) + 2N < +cx3. 

We claim that for any r S z^^, we also have 

(4.29) JUXe„,Q1) > JUXe„^^ + QI) - en 

for some error e„ — > which is independent on t. 

Assume to simplify that r = svi for some s e M. If s is (1 — 277) times an integer, then (|4.29p 
is obvious. If not, we may assume without loss of generality that < s < 1 — 277. We have 

JhMe ^s^i + QD = T- F{xe [x + hnJ:])dx. 

The domain of integration is spht into three parts Di = l''" {svi +Q2)n (Q^), -D2 — I^" {svi + 
Ql) n /''"((I - 277)z.i + Ql), and = I^^-isi^i + Ql) \ {D, U D2). Since Di and {D2 - (1 - 277)1^1) 
are disjoint subsets of I^"{Q2)- it follows 

(4.30) JhAxE.^su,+Q2) < JUxe.^QI) + ^ f F{xE\x + hnn)dx. 




JD3 

We have £13 C 5"!+ U ((1 - 277)1/1 + S:[), where for t = 1, . . . , iV - 1, 

N 

Let us show that (l//i„) Jg+ F(xe [^ + Ki^]) dx —^ as n —> 00. We have, using the notation 
X = siyi + x' , x' ■ vi ^ 0, and the change of variable s = 1/2 — 77 — /i„^, 

(4.31) / F{xe [x + hnn)dx = / F{vn{£.,x'))dx' d^ 

iT-n JS+ " JQ J[-l/2+i),l/2-i)]»-i 

where for each ^, x', 7;„(^, x') G {0, 1}^ is the vector [(1/2 — 77— /inO^i+^' + ^'^n^]- We observe 
that from (|4.26p . we have (using the same change of variable, and observing that x/„ depends only 
on x'jv) 

/ / \xiAx'^)-XEA^/2-V-hni,x)\dx'd( ^ 

Jo J[-l/2,l/2]«-i 

as 71 — > 00. In particular, up to a subsequence, we may assume that each component of the 
vector Vn{£.,x') converges to Xi^ix^) as n — > 00, for a.e. (^,a:'). Notice that these components 
take only the value or 1. hence, they must be equal to X/,^(^Ar) ^'^^ large enough. Since 
F{0) = F(xe) = 0, it follows that F{vn{£,,x')) a.e., and since F is bounded we find (using 
Lebesgue's theorem) that the integral in the right-hand side of (|4.31|) goes to 0. 
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In the same way, wc can show that (l//i„) Jg- F{xj^ [x + /in^]) da; 0, so that (|4.29p foUows 
from (|4.30p and the inclusion 1)3 C Si U ((1 — 2ri)vi + S^). In the general case (if r is not parallel 
to ui), wc can show in the same way that ()4.29p holds provided 

our choice of 77 guarantees that e„ still goes to zero. 

Let Un '■ Q2 ~^ [Oj 1] be the average of on each hyperplane orthogonal to j^, given by 

It is clear that w„ depends only on a; • v. Since 

/ \un- XiAdx < / \xe -XiJdx, 

we deduce 

Notice that for a.e. a; G Q2, the vector Un[x + is also the average on the hyperplane through 
X orthogonal to of the vectors Xe [' + so that, by Jensen's inequality and using (|4.29p 

Together with (|07)) . it yields 
(4.32) e > liminfj;; iun,Q1). 

n — yoQ " 

It is clear that (|4.28p also yields a uniform bound on the total variations \Dun\{Q2), n G N. 
Let us fix e,5 e]0, 1[. By the generalized coarea formula for , and the coarea formula in BV 
we get 

C > JU^n,Ql)+e\Dun\{Q2) = f\j;iJx{u^>s},Q2)+en''-\d*{un > s}nQl)]ds 







> 



for some positive constant C. We deduce that there exists s„ g]S,1 — 6[ such that 

(4.33) JfiJun,Q2)+e\Dun\iQ'l) > (1 - 2^) J,^Jx{n„>s„}, Q^) 
and 

(4.34) C > £(1 - 2S)n''-\d*{un > s„} n Q]>). 
We have 

Vn ■■= X{u„>s„} Xiu in L'^iQ'D- 
By (|4.34p and since w„ depends only on a; • we deduce that there exists, up to a subsequence, 
an odd number M independent of n such that 

M 

5*K = 1} = QI! n |J(af;;. + i^^) 

with 

-^+V < al < al < ... < < ^-V- 

Without loss of generality, we can assume that for every A: = 1,...,M, w„=Oifa;-j/<a^ 

and w„ = 1 if X • > . Indeed, if some of the points do not go to zero as n — *■ 00, we can 
lower the energy (wn,QU) by "removing" from Vn the corresponding discontinuities. 
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If M > 1, let US consider the function Vn ~ X{x iy>al}- Since w„ V w„ is both a translate of Vn 
(of (a^ — a^)i^) and of xi„ (of o.n'^), for n large we have 

By the submodularity property of F, we obtain setting := u„ A u„ that for n large enough, 

If A/ = 3 then = X{x i^>a->} is a translate of xi„- If > 3, then we can reiterate this argument, 
replacing now with a^: after a finite number of steps we eventually find a translate of xi^, with 
energy lower than J^^{vn,Q2)- Hence, for n large enough, 

(4.35) JUxi^^QD < JLi^n,Q1). 

By (|4.16p and a straightforward scaling argument, the left-hand side of (j4.35p goes to (1 — 
2rf)^-'^F{v • E) as n ^ oo. Recalling (j4.32p and (|4.33p . we deduce 

l + esn^\Dun\{Ql) > (1 - 2S){1 - 2i^f-' F{iy ^) , 

and since e, 6 and rj can be chosen arbitrarily small we deduce that £ > F{v ■ S). It follows that 
fi>F{i^- T.)H^-'^Ld*E, which together with ([4T9ll yields (j4T^ . 

4.2. F-limsup inequality. Let us now show the inequality (|4.13p . We first show the following 
generalization of the formula (|4.16p : 

Lemma 4.6. Let A C be an open set, and assume ^.^^^{dly n dA) — 0. Then 

(4.36) lim j,^(x/„,^) = n''-\dhr\A)F{u-i:). 

Proof. Let e > 0. Let (j^Oi^i be a basis of K^, with = v, and let us consider the family Qe of 
all cubes Q C ^, of side e, centered at the points X^ili ^^i^i for fc^ G Z^, and such that each face 
is orthogonal to some Vi. 
We have obviously 

lim J^(X/,.,^) > Yl liminf J^(x/„,Q) 
while (by (j4.16p and a simple scaling argument) for each Q, 



Hence, 



lim4^(x/„,g) = e''-^F{vY.). 



lim J^h{xi,.,A) > n^-' idl^n y Q]Fiiy E) 

QSQa: 



letting then e ^ 0, it shows ">" in (|436l) . 

To show the reverse inequality, we let, for each e > 0, = Uqgq ^.nd first observe that 
for any 77 > 0, 

Jh{xi.:A-) < J,^(x/„,(l + r?)g) 

as soon as ft. < rie/{2p-s), where here (1 + t])Q denotes the cube of same center as Q, and dilated 
by the factor 1 + 77. Taking the limit ft ^ 0, we find 

(4.37) limsupJ^(x/,,A'^) < il+v f-'^n'^'^idIr,nA')F{iyj:) < l-L^~\dI^r\A)F{u-Y.) + Cr^ . 
for any t] > 0. Now, 

Jh{xi..:A) - r^{xi.,A') + \ [ F{xiAx + hT])dx 
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SO that (|4.36p will follow from (|4.37p if wc show that 
(4.38) i / F{xiAx + hJ:])dx ^ 

" ■'/''(A)\/''(A=) 



as /i — > and e — > 0. The integrand above is zero when x is at distance larger than hp^: to the 
interface dli,, and x € A is out of the domain of integration as soon as it is at distance larger than 
Zips to A \ , for instance when x G := G ^ ■ dist(^, 9^) > 2e}. Hence, the error (|4.38p 
is bounded by 

(4.39) (maxF)'^-^\-^^-^^^^^- + [-"^^'^^^]'-)'. 

E ' h 

Let G C dli, be a relatively open set which contains dli, D dA. We claim that there exists (5o > 
small such that the projection of {A \ As) n (9/^ + [—5,6]i^) onto dl^, is contained in G for any 
S < So', if not, one finds a sequence (S„ and points Xn € {A \ Ag^) O {dli, + [—Sn,Sn]v) which 
project outside of G, but then, any limit point of this sequence should be in dA n dli, but outside 
G, which is not possible. 

Then, if we choose e < So/2\^, and h small enough, we have 

{A\A2e)n{dI^ + [-hp^,hp^]iy) c G+[-hpj:,hp^]u 

so that the limsup of (|4.39p . as /i ^ 0, is less than 

PS (maxF) H^-\G). 

Since we assumed that 'H^~^{dl^ n dA) = 0, Ti.^~^{G) may be chosen as small as we wish so 
that (|4.38p holds. Hence the Lemma is proved. □ 

Now, we show the following estimate: 

Lemma 4.7. Let A d fl be an open set with continuous boundary, and let E d fl be a finite- 
perimeter set such that n^~'^{d*EnA \ d*E n A) = 0. Then, 

(4.40) limsup J^(x£;, A) < {2psmaxF)n^^\d* E n A) . 

Proof. Since A has a continuous boundary, for h small enough, x + hYi contains points of Lcbesgue 
density of E both zero and one only if x is close enough to d*E{^A, namely, dist [x, d*E{^A) < hp^. 
Hence, 

jc( ,^ 1 /• p. r ^-viA^ <- : dist{x,d*EnA)<hps}\ ^ 
JhiXE.A) = - / F{xe[x + hi:\)dx < (maxF). 

By standard results on the Minkowski contents [9^ Thm 3.2.39], the last fraction goes to 
2pYH^^'^{^*Er^ A) as /i ^ 0, which shows the Lemma. □ 

We deduce the following: 

Corollary 4.8. Let Vl C have a continuous boundary, and let E (Z fl be a set whose boundary 
is made of a finite union of subsets of (N — I) -dimensional hyperplanes. Then 

\imJf;{xE,n) = f F{vE{x)-^)dn''-\x), 
''^0 JdEnn 

where i^eIx) is the inner normal to E at x. 

Proof. By assumption, dE nVl — Uf=i where Pi C {xi + 9/^.) for some Xi e M.^ and e 
with unit norm, moreover we assume that i^i = ve on Pi {vi points towards the interior of E). (We 
also assume that the Pi are "maximal", in the sense that Pi Pj ^ 9 ^ Vi ^ Vj for any i ^ j.) 
Let 5 = n Ui^f^i ^^i^ where the dPi is the relative boundary of the face inside the hypcrplane 
{xi + dl^-): it is a (iV — 2)-dimensional set, with H^~'^{S) < +00. 
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Wc choose £ > and build a finite coverine; 

iA)iit^ of dE nn with bounded open sets 
with continuous boundary, as follows: = {x G : dist (x, i9r2) < e}, Am+2 = {x £ f^, : 

dist (x, S) < e}, and Ai is a neighborhood of Pi \ {Anj^^i U Anj^2) which does not intersect Pj for 
j i- By Lemmas 14.61 and 14. 7[ 

AI 



lim r^ixE,^) - y2n^-'{P; n A,)F{iy, ■ E) 



< c {n''-\dE n Am+i) + h^^^oe n Am+2)) 



with C ~ 2pE maxs F. Letting e ^ shows the Corollary. 



□ 



We are now able to show the following Proposition, which essentially shows the F-convergence 
of J/, to J on the restricted class of polyhedral sets. 

Proposition 4.9. Let il C have a continuous boundary, and let E (Z ^ be a set whose 
boundary is made of a finite union of subsets of (N — 1) -dimensional hyperplanes. Then, there 
exist sets Eh with xEh ^ Vh{^), XEh ^ Xe "in L^{^) as h 0, and 



liia JhixE,,,^) 



dEnn 



F{Mx)-^)dH^-\x) = J{xE,n)^ 



Proof. We have, making in (|4.14p the change of variable x = z + y with z G /iZ^ and y G [0,h[^ , 



Jh{XE,n) 



1 



[OM 



dy 



h""-' E F{xE[y + z + hn 

so that (by Corollary 14. 8p we can choose for each h a, yt & [0, h[^ such that 
(4.41) hmsup /i^-i F{xE[yh + z + hn) < J{xE,n), 

we can assume moreover that no point in yh + hl^^ H il lies on dE. 

For each e > we let fi^ = {x e il : dist (x, d^l) > e}. We fix e > 0, and define Ef^ as follows 
(for h small enough): 

XB^ = Y XE{yh + z)xQii + Y Xe{z)xqi^ 

where Q'^ is the cube defined in (|4.3p . It is not difficult to show that XEf^ Xe in as h ^ 0, 
in fact, it converges locally uniformly in \ dE. We have if h is small enough 



(4.42) 

A{xEi,n) < h^-' Y F{xEi[x + hj:]) + h^-^ Y F{xEf}x + hi:]) 
</i^-i Y F{xE[yh + z + hn) + h''-^ Y FixE'Jx + hi:]) 

ze{i'^{n)~yh)nhz'^ a:e/''(a\f22e)n/iZ" 
Now, thanks to (|4.41[) we get 

(4.43) limsnpJh{xEf^,n)<JixE,n)+limsuph^'^ ^ F(x£j [x + /lE]). 

xei'^{n\n2,)nh'L« 

In the sum, on the other hand, F{xei [x + hT,]) is not zero only when some point of x + hY, lies in 
Ef^ and some other in il\Ef^, and such x are at distance at most /i(l + ps) from dE n il, so that 

c f^25+ftV]v • dist (e, dEnn)<h{i + ^ + p^)} 

Therefore, the sum is bounded by 

lUe »\»2e+/.viv • distic,dEnn)<h{i + VN + p^)}\ 

h 
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which goes to 2(1 + \/N + ps)H'^~^{dE D {n\ fJaJ) as /i ^ 0. Wc deduce from that 

Since TC^~^{dE \ il2e)) — > as e — > 0, using a diagonal argument, we deduce the thesis of 
the Corollary. □ 

Corollary 4.10. Let SI C M.^ have a continuous boundary, and let E be a polyhedral set in SI in 
the previous sense. Then, 

(r- lim A(.,r!))(xB) - J{xE,n). 

Proof. It follows from (|4.12[) (which has been shown in Subsection l4.1l) and from Proposition 14.91 

□ 

Now, we are in a position to show that ()4.13p holds. 

Proposition 4.11. Let il be a bounded open set ofR^ with Lipschitz boundary, and let E C ft 
be a set with finite perimeter in fl. Then for every h > there exists Eh with xe^ G Vh{^l), such 
that XEh ~* Xe in L^{Q,) as h and 

(4.44) limsupJ,,(xi5,,S7) < J{xeM)- 

Proof. Since S7 is Lipschitz, we can extend E outside of f2 into a finite-perimeter subset of M.^ 
(still denoted E) such that \DxE\{di^) = 0. Then, standard approximation arguments show that 
there exists a sequence of polyhedral sets G" such that xg" converges to xe strongly in L^(]R^), 
and with 

hm \DxgA^'')\ = iDxEm""). 

n — >oc 

This can be seen, for instance, by approximating xe by smooth functions (by convolution) and 
then approximating these smooth functions by piecewise linear functions, such as "PI" finite- 
elements. Then, an appropriate thresholding of these functions provides the sequence G". The 
Reshetnyak continuity Theorem (see section[2|), and \DxE\{d^) — yield 

lim J{xG",^) = J{xE,^)- 

n — >oo 

By Proposition 14. 9i we know that for each n there exists Gf^ converging to G", such that 
limsup^^Q J?i(GJJ, ri) < J{xg"t^)- We construct the family {Eh)h>o from the G^, by a diag- 
onal argument as follows. For every n there exists such that /i„ | as n ^ cx3 such that for 
every h < hn we have ||xg;; - XG"||Li(n) < 1/"- and JhixGi,^) < JiXG^,^) + ^jn. If we set 
Eh G\ for < h < hn, the result follows. □ 

5. Examples 

Let us describe a few cases. First of all, the standard nearest-neighbor interaction on a square 
grid corresponds to the situation where S = {(0, 0), (1, 0), (0, 1)} and, for u G M^, F{u) = \u{l, 0) — 
m(0, 0)1 -1- \u{Q, 1) — m(0, 0)|. It is obvious, in this case (as in any other case where is a sum of 
pair interactions) that the F- limit of Jh is the anisotropic total variation given by (|1.5p . in this 
case, \Du\i where \p\i = \pi \ + \p2\ is the 1-norm in M^. 

Less trivial situations are when F cannot be reduced to a sum of pair interactions, such as, still 
with the same set E, the functions F defined by i^(Os) = -F(ls) = 0, F{1^ — u) = F{u) for any 
u € {0, 1}^, and 

f(' Ul, f(' Ul, f(' ]=V2. 

\ 1 ' voo; ' vol' 



22 



A. CHAMBOLLE, A. GIACOMINL AND L. LUSSARDI 



This F can also be checked to be submodular. Now. the limit density is given by 

1^2 



(5.1) 



F 



vi 

see Figure[T]for the expression of F, and where we also have plotted the shape of the "Wulff shape" 
F{p • E) < 1. Notice that in this case, we have chosen F{9) = ^/(f^i.o — ^o,o)^ + (6*0,1 — Oa.o)^ for 



V2\V2\ + 



V2\U2 \ + W2 




+ — 



/ V2Wl\ + W2 - I'll 



Figure 1 . Values of F given by (|5.ip 

9 e 1]^°'^^: however, clearly, we get a limit energy which is not /^^ \Du\ (which would be the F-limit 
of u /i - \/(wi+i.j — "^ijY + ('^i.j+i ^ "^ijY) but an anisotropic (crystalline) total variation. 
This is of course due to the fact that the latter does not satisfy the generalized coarea formula. 
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